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Abstract 

We model bond’s price curves corresponding to the sovereign Uru¬ 
guayan debt nominated in USD, as an alternative to the official bond 
prices publication released by the Central Bank of Uruguay (CBU). 
Four different gaussian models are fitted, based on historical data is¬ 
sued by the CBU, corresponding to some of the more frequently traded 
bonds. The main difficulty we approach is the absence of liquidity in 
the bond market. Nevertheless the adjustment is relatively good, giv¬ 
ing the possibility of non-arbitrage pricing of the whole family of non 
traded instruments, and also the possibility of pricing derivative secu¬ 
rities. 


1 Introduction 

Bond prices curves (or equivalently yield curves) constitute a major tool 
in debt analysis and perspective of the sovereign debt of a country. Term- 
structure models have therefore been used in different ways by different 
classes of market participants. In the monetary policy context, the term 
structure is an indicator of the market’s expectations regarding interest rate 
and inflation rates. From a financial point of view, the existence of a bond 
price curve helps the development of the domestic capital market, both for 
primary and secondary market. 

There are essentially two approaches to model of the term structure. The 
general equilibrium approach starts from a description of the economy and 
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derives the term structure of interest rate endogenously, as used for example 
by Cox Ingersoll and Ross in [0j. In contrast the arbitrage approach starts 
from assumptions about the stochastic evolution of one or more interest 
rates and derives prices of contingent claims by imposing the no arbitrage 
condition, this is, for example, the pioneering approach proposed by Vasicek 

Ca¬ 
in this first approximation to the problem of the Uruguayan USD nomi¬ 
nated debt, we use the second approach and further we restrict our analysis 
to four gaussian models because of their analytical and numerical tractabil- 
ity. Factor models assume that the term structure of interest rate is driven 
by a set of “state variables” named as “factors”. As one factor generally 
explains a large proportion of the yield curve movement, it may tempting 
to reduce the analysis to one factor models. Nevertheless, the consideration 
of one factor implies the correlation between two different rates in the same 
time interval. This perfect correlation is difficult to accept as empirical data 
usually shows high correlation, so then we use multiple factors model, and 
in particular, two factor models. 

The first two models we choose are short rate models, the classical Va¬ 
sicek model [13] . and the more flexible G2+-1- model [2]. The parameter 
estimation is carried out based on maximum likelihood, following Chen and 
Scott [3j. The second couple of models belong to the HJM family [Tj. Here 
we again choose first the simplest possible one, the Ho-Lee model with con¬ 
stant volatility [8], and second, the more flexible Hull-White model with 
tempered volatility [5]. In this two last cases we calibrate the models with 
the help of minimization of squared differences. 

Our contribution is in first place to provide an arbitrage-free set of bond 
prices for the Uruguayan USD debt, that can be used to portfolio valuation 
and derivative pricing. 

The rest of the paper is organized as follows. Section 2 provides infor¬ 
mation about the Uruguayan debt and the publication of prices by CBU. 
Section 3 describes the short rate model we use in the work and and the 
estimation methodology through maximum likelihood. In Section 4 we de¬ 
scribe the results in the Vasicek and G2+-1- models. In Section 5 the HJM 
models are introduced, and the calibration procedure is explained. Section 
6 describes the results in Ho-Lee and Hull-White models, and Section 7 
concludes. 
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2 Bond price curves 

2.1 About the Uruguayan debt 

The Uruguayan government issues debt through different financial instru¬ 
ments, with a variation of currencies and expirations. In present times, the 
most relevant circulating instruments according to terms and currencies are 


• the Treasury Notes in local currency and linked to CPI, 

• the Local Bonds (issued in the country) linked to CPI (consumer price 
index) and in USD, 

• the Global Bonds (issued mainly in United States of America) in USD. 


To complete the analysis of the debt, it should be mentioned that the 
government also holds obligations with multilateral financial institutions, 
estimated in less than the 10% of the global debt amount. In Figure 2.1 
we observe the debt’s profile of the Uruguayan debt across maturities. The 
last expiration dates correspond to the 2050 bond, recently issued, paying 
its face value on thirds, in the years 2048, 2049 and 2050. 
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Figure 1: Central Government Debt Profile 

Regarding the currency composition, more than one half of the Uruguayan 
debt is nominated in local currency (part of it linked to CPI), 45% is issued 
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in USD, and the rest corresponds to bonds issued in Euros and Yens. Al¬ 
most all of the debt (94 %) pays fixed coupons, the rest is mainly adjusted 
to Libor rate. 

2.2 Curve price estimation 

In most of the countries the corresponding regulating agencies and some 
financial corporations release yield curves corresponding to the sovereign 
debt. The most popular methods include parametric methods (as the ones 
proposed by Nelson-Siegel HU and Svensson |12| ) also used non-parametric 
methods lesser extent. 

Regarding Uruguay, the CBU does not release yield curves, but daily 
issues a report with prices for all the financial instruments issued by the 
government. The two authorized stock exchange institutions release their 
respective yield curves. The Electronic Stock Exchange of Uruguay (Bolsa 
Electronica de Valores del Uruguay - BEVSA) uses B-splines to produce the 
yield curves corresponding to local currency, linked to CPI, and to USD. 
The Montevideo Stock Exchange (Bolsa de Valores de Montevideo - BVM) 
emits a single curve (linked to CPI) based on the methodology proposed by 
Svensson jl2 j. 

The purpose of the present paper is to propose an alternative methodol¬ 
ogy for bond valuation, for the debt issued in USD, to the one used by the 
CBU and also by the BVM and BEVSA. In the next subsection we briefly 
describe the procedure used by the monetary authority (CBU) to compute 
the prices of the different financial instruments. 

2.3 CBU pricing methodology 

The CBU issues daily reports containing prices of all the financial instru¬ 
ments issued by the Uruguayan government, named as “price vector” and 
available in the local stock exchanges. The pension funds and insurance 
companies, that are relevant participants in the domestic bond market, are 
obliged by law to use this prices in order to report portfolio values. 

This vector of prices is computed through a methodology that includes 
four different criteria, according to whether the bond has been traded or not, 
and according also to its expiration, and apply with a hierarchical scheme. 
At the end of each business day the CBU releases information about bond 
prices, according to the following procedure: 

1. For bonds that have been negotiated in the day (according to certain 
minimal amounts) the prices are computed as a weighted mean of the 
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respective negotiation prices in the stocks exchanges. 

2. For bonds with less than a year to expire, an interpolation procedure 
between the previous price and the face value is applied to compute 
the prices 

3. An index I(t ) where t represent the current date is calculated in order 
to compute the prices of the other bonds. This is a number that 
represent the value of an “ideal” mean bond. 

4. For all other bonds, the new price is the previous price multiplied by 

the ratio computed. 

The most relevant characteristic to be taken into account when analysing 
this procedure, is that the stock exchange of long expiration debt instru¬ 
ments issued by the government (bonds) is really not liquid. This implies 
that most prices published by the CBU are computed instead of negoti¬ 
ated (see points 3 and 4 above), and this can happen for some instruments 
consecutively, during a relatively long period of time. This can lead to (the¬ 
oretical) arbitrage opportunities, for instance giving larger yields for bonds 
with smaller maturities, what is equivalent to larger prices of zero coupon 
bond with smaller maturity than other zero coupon bonds. Of course, this 
arbitrage is eliminated in real negotiation, but is present, for instance, when 
evaluating portfolios. Although the intention of the methodology is to follow 
the movements of the market, one should always take into account this fact. 

3 Short rate models 

When modelling yield curves through stochastic processes, the classical ap¬ 
proach consists in modelling the short rate through a certain amount of 
sources of uncertainty, under the denomination of “states”. When this states 
are used to construct the short interest rate we obtain a large variety of mod¬ 
els, described in a vast part of the literature in fixed income mathematical 
finance m, m, u and ma¬ 
in this part of the work we use two gaussian models, the one proposed by 
Vasicek in 1977 E3I, that is seminal in the literature. The second model we 
apply is the G2++ [2], that is a modification the previous one that has two 
factors and includes the initial price curve, avoiding in this way arbitrage 
at the initial time. This model has in fact a more general version, as it can 
include an arbitrary amount of factors, and is called the Gn+-1- model [5] 
in this case. However, the choice of the number of factors then involves a 
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compromise between numerically efficient implementation and capability of 
the model to represent realistic correlation patterns and to fit satisfactorily 
enough market data. 

In the next two subsection we review each of the models to be used. 

3.1 Vasicek’s model 

In his classical paper [13], Vasicek proposes a model for the short rate 
through a stochastic differential equation driven by a Wiener process, 

dr(t) = a(b — r{t))dt + crdW(t) r( 0) = ro; 

where a, by a are positive constants and W ( t ) is a standard Wiener process. 
The solution to this equation is know as the Ornstein-Uhlenbeck process. 
It defines an elastic random walk around a trend, with a mean reverting 
characteristic. Given the set of information at time s , the short rate r(t ) is 
normally distributed with 

E (r t \F s ) =r s e ~ a + b(l - e 

Var(r t \F s ) ~(l -e" 2 ^)). 

The bond price can be obtained by computing the discounted expected 
terminal value of the bond with respect to a risk neutral probability measure 
Q. This quantity can be explicitly computed concluding that the Vasicek 
model is an affine model whose the solution is 

P{t, T) = A(t, T)e ~ B{t ’ T)n ; (1) 


where 

2 2 

A(t, T) = exp ((6 - 2 ) (B(t, T)-T + t)~ £ B(t , T) 2 ); 

B(t, T) = 

3.2 The G2H—b model 

As mentioned above, the price correlation given by the Vasicek model to 
different instruments is 1. For this reason, as observed prices do not show 
such high correlation, emerged the idea to introduce models with more fac¬ 
tors, that imply more parameters, that would allow to better fit to the data. 
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One of this proposals is the Gn++ model, proposed in [5] that we use with 
two factors. Another relevant characteristic of this proposal is that in takes 
into account the whole initial price curve. 

The dynamics of the short rate process in this model is given by 

r(t ) = x(t) + y(t) + ip(t), r(0) = r 0 ; 

where the process x(t) y y(t) and the vector (ITi, W 2 ) satisfy 

dx = —axdt + adW \; 
dy = -bydt + rjdW 2 ] 
dW\ (t)dW 2 (t) = pdt ; 

where ro,a,b,a,r] are positive constant and — 1 < p < 1. Given the set of 
information at time s, the short rate r(t) is normally distributed with 

E(r{t)\F S ) =x(s)e~ a ^ + y(s)e ~ b <*-> + ^(t); 

Var(r(t)\T s )=^(l-e~ 2a ^) 

2a 

+ |!(i _ e -m-s)) + ^^(1 - ex P -( a+b ^- s )). 

The price at time t of a zero coupon bond with maturity at time T is 


P(t, T ) = exp ^ ^ <p(u)d 


1 _ e -a(T-t) 

u --- x(t) 

a 

1 _ e -b(T-t) 


y(t)+ l -V(t,T))- (2) 


where 


a 


V (t, T) = — (T — t -\—e 


2 0 -a(T-t) _ J_ e -2 a(T-t) 


2 a 


3 

2 a 


+ r L(r-t + 2 e - fe (' r - t ) - Jig-2 b{T-t) _ 

+ b 2 V t+ b 2b 2b) 


2parf_/ 
ab V 


-a(T-t) _ ^ -b(T-t ) _ y e -(o+6)(T-t) _ ^ 

T - t +-+---- 


a 


a + b 


(3) 


Even though the previous formula ([ 2 ]) gives bond prices in the model, it 
is necessary to estimate the function p>. In order to do this, it is necessary to 
assume that the initial price curve T 1 —> P AI (0,T) is known. The model fits 
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the currently observed term structure if p Model ( o, T) = P M ( 0, T), therefore 
the price at time t of a zero coupon bond maturity at time T is 


P(t,T) 


P m (0,T ) 
P M (0, t) 


exp(A(t,T)); 


(4) 


where 

1 \ i _ e -a(T-t) i _ e ~ b ( T ~t) 

A{t, T) = -( V (t,T)-V(0,T)+V(0,t)) --- x(t) --- y(t). 

3.3 Maximum Likelihood Estimation 

Observe that in formulas 0 and the log-prices are expressed as a linear 
function of the state variables. In order to determine this state variables in 
both models it is necessary to use one bond price time series in Vasicek model 
and two time series in G2+-K But to carry our this procedure the values 
of the parameters (to be estimated) is necessary. This is why we use the 
Maximum Likelihood (ML) method with the help of a change of variables, 
that gives the corresponding Jacobian term in the ML expression. In both 
models, based on the Markov property of the processes, the joint density 
is written as a product of conditional densities, each of one has normal 
distribution, with parameters three parameters in the first case (a, b, a), 
and five in the second (a, b, a, 77, p). 

Given the panel data set P t M = (P M (i, P*^)), i = 1y t = 1,...,T, 
where P M (t,T^) is the price at time t of the zero coupon bond with ma¬ 
turity pW and denote by X t = (x^), t = 1,..., T the state vector. 

The joint density of P ^ 1 ,..., P^ ! satisfies 

/(P 2 m ,...,P^|Pi M , 0) = n/(^l^ 0 ) 

k =2 

by the Markov property. Changing variables in each conditional density, we 
obtain 

/(Pf iP^r,©) = f(X k 0)1^1 

= /(X fe |X fc _ 1 ,0)|^0| = /(X fc |X fc _i,©)^; 

dx k 1 fc| 

where Jf. is the Jacobian of the change of variables. 
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Remark. We find adequate in our situation to use as many bond price time 
series as number of factors. But this is not strictly necessary. In some cases, 
when more time series than state variables are used it is possible to introduce 
additional random variables in order to estimate the model [3]. 

4 Short rate estimation 

4.1 Data 

The information used in the construction of the yield curves is the one 
provided by the CBU, for USD nominated bond prices, traded both in the 
domestic market an in foreign exchanges. It should be noticed that, besides 
the liquidity problem of the bond market, there are no derivatives on this 
instruments, so all the available information is provided by the bond prices. 

The data to be used in the estimation procedure corresponds to the time 
period January 4, 2010 to October 30, 2013. In the Vasicek model we use 
the bond BE330115P, expiring in January 2033, and in the G2++ model 
we use the same one an add the bond BE250928F with expiry in September 
2025. We choose these bonds as they are the more frequently traded. Each 
time series is processed according the coupon payment scheme (amount and 
frequency) to obtain the corresponding yields, that give the zero coupon 
bond prices used in the estimation. 

4.2 Estimation in the Vasicek model 

In order to evaluate the influence of the non liquid computed prices (as 
explained above) we proceed in this case estimating the parameters corre¬ 
sponding to the Vasicek model in two situations. We first estimate parame¬ 
ters using all the available information, i.e. we use weekly bond prices taken 
on Wednesdays. The obtained parameters are: 

o = 1.7051; 6 = 0.0937; a = 0,3721. 

In the second case we use only negotiated prices (504 observations). This 
implies that the time intervals between prices are not regular, leading to a 
slightly more complicated estimation scheme. The obtained parameters are 

0 = 1,7145; 6 = 0,0896; a = 0,4971. 

This comparison help us to verify that the variation in the values of the 
parameters is not significant (perhaps the most important variation is reg¬ 
istered in a) and for this reason in what follows, and in this first approach 
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to the problem, we will use all the available time series. In Figure 4.2 


we 


observe the price curve corresponding to August 13, 2013 given by the model 
in both estimations carried out, and also the bond prices issued by the CBU. 



Figure 2: The solid line shows the first estimation with Vasicek model, 
the dotted one the second (that uses only negotiated prices). The dots 
correspond to market prices of zero coupon bonds. 


4.3 Estimation in the G2H—b model 

As we mentioned above we use two bond prices time series, with expira¬ 
tion in 2025 and 2033 respectively. We take weekly prices corresponding to 
Wednesdays. The obtained values with the ML estimation described above 
are 


a = 0.1300; b = 0.3526; a = 0.2062; 77 = 0.4892; p = -0.99. 

Using this values with the corresponding bond prices we can obtain the daily 
time series corresponding to the short rate. This allow us to compute the 
bond prices for arbitrary maturities in each one of the days. We show the 
bond price curve for the following maturities: 1, 2, 3, 6 and 9 months and 
1, 2, 3, 5, 7, 10, 15, 20 and 25 years. 

In Figure [T3| we see the daily zero coupon bond price curve corresponding 
to the analysed time period issued by the Uruguayan government, as a result 
of the application of the G2+-(- model. 

If we analyse the data in Figure |4.3| for the second semester of 2012, 
the model gives a curve with non-negative slope for some maturities. This 
fact introduces arbitrage possibilities, as it gives cheaper bonds with smaller 
maturities than others with larger maturities. More precisely, our model 
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Figure 3: Surface of zero coupon bond prices using the G2++ model 


adjustment give parameter that result in increasing bond prices in some 
intervals of the maturity T. Our parameter set includes a correlation very 
close to —1, and instances of negative values of the factors x(t) and y(t). 
We will return to this fact in more detail in the next subsection. 

4.4 On arbitrage possibilities in the G2H—|- model 

We analyse how formula Q for a zero coupon bond price in G2++ model 
can give arbitrage opportunities. Equation (|3]) can be written as 

V(t, T) = ^l [ (1 - e~ a( ' T ~ u ' ) ) 2 du + [ (1 - e- b(T ~ u) ) 2 du 

a 2 J t b 2 J t 

+ 2pm 1 r T (1 _ e _ a{T _ u) _ b{T _ u) 

ab J t 
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We then take logarithms and differentiate formula Q w.r.t T, to obtian 
^P(t, T) = ^P M ( 0, T) + ^ ((1 - e-“M) 2 - (1 - e-“ T ) 2 ) 

+ ^( (1 - e " 6(T " t))2 - (1 - e " 6T)2 ) 

+ ^ ((1 - e-“M)(l - e _fe ( T_f )) _ (i _ e -“ r )(1 - e - feT )) 


In this formula the first three addends are negative. The fourth has the op¬ 
posite sign of p, the last two depend on x(t) and y(t). In our application p is 
close to —1, what is associated with the fact that x(t) and y(t) take opposite 
signs, therefore, in some times intervals the derivative is positive, as shown 
in Figure 4.3 Nevertheless the model gives general valuable information 


about the structure of the debt. 


5 Heath-Jarrow-Morton model 

In the market we do not have a real instantaneous interest rate. In certain 
cases, the one (or three) month interest rate series is used as a proxy to 
estimate this interest rate. It is not convenient to use the overnight rate, as 
it has very high volatility due to economical factors, as for instance liquidity, 
and this can not depend on the structure of the yield curve. 

By this reason, and also to avoid arbitrage opportunities in a systematic 
way, Heath, Jarrow and Morton introduce a new methodology [7] (refereed 
as HJM models), that models the forward instantaneous rate at time t by a 
stochastic differential equation driven by a Wiener process 

f(s,t) = f(0,t) + f a(u, t)du + [ a(u, t)dW u ; 

Jo Jo 

where Wt is a Wiener process. 

We depart from the free arbitrage model, that assumes that the a(u,t), 
0 < u < t and a(u,t), 0 < u < t are adapted processes defined in an 
underlying stochastic basis (Q, J 7 , P u , Q ), and Q is a risk-neutral probability 
measure. The key aspect of HJM techniques lies in the recognition that the 
drifts of the no-arbitrage evolution of certain variables can be expressed as 
functions of their volatilities and the correlations among themselves. 

a(s,t) = <r(s,t) / a(s,u)du. (5) 

J S 
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Therefore, in order to specify an HJM model the initial forward rate 
curve f(0,t) and the volatility structure a(s,t) should be given because 
no drift estimation is needed. It should be observed that as long as the 
function a is deterministic, by condition © the drift is also deterministic, 
in consequence, the forward rates are gaussian. 

In the applications that follow, we first choose the simplest possible 
alternative for the volatility, i.e. we assume a(s,u) = o > 0 is a positive 
constant. This gives 

f(s, t) = /(0, t) + a(st - s 2 / 2) + aWt) 

the so called Ho-Lee model |8j. The price of zero coupon bond in this model 
is then given by 

P(t , T) = P pJ ( ° Q J ) ) exp ((T - t)f M { 0, t) - y t(T - t ) 2 - (T - t)r(tj) . (6) 

In the second model we assume that cr(s,u) = ae~ a ^ u ~ s \ where a and a 
are positive constants. This gives 

r(t) = /(0, t) + - l) 2 + J* e-^-^dWu. 

The price of a zero coupon bond in this case is 

P(t, T ) = P pJ^ exp (5(t, T)/ M ( 0, t) 

~ ^(1 - e~ 2t )B(t, T ) 2 - B(t, T)r{t )). (7) 

that corresponds to the Hull-White model [9] with time dependent param¬ 
eters. 

6 Results for HJM models 

In order to use it in practice to price contingent claims, we have to calibrate 
its parameters on market data. To do this we can adopt a cross-sectional 
approach. Assume we observe a cross-section of market prices of contingent 
claims, that is to say, the prices of a set of N zero coupon bonds made by 
the market at the same day. For time t, let us denote the vector prices 
P t M = i = 1 ,...,/ where P M (t,T^) is the price in time t 
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the zero coupon bond maturity at time T ® and assume we are able to 
compute the price vector made by the model and denote them by p^ odel = 
(pModel^ , i = 1Model prices are functions of the parameter 

vector of the model, see © in Ho-Lee and ([7]) in Hull-White. 

The idea is to find the values for 0 minimizing the difference between 
market prices and model prices. To do this we have to solve the least square 
problem, 

min V (P M (f,T W ) - pModel^ T (i)^y 

i— 1 

In this analysis we use weekly prices corresponding to the whole year 
2014, of 10 USD nominated bonds issued by the CBU. The maximum ex¬ 
piration date corresponds to the BE451120F bond, in November 2045. To 
obtain an approximation of the initial interest rate we used the reference 
curve CUD-BEVSA daily issued by BEVSA for three months, and to ad¬ 
just the initial forward rate curve we used the yield curve corresponding to 
January 5, 2013. 

Calibrating the Ho-Lee model we obtain a a for each day. In Figure |4] 
we observe then the daily variation of a. The mean value of the estimation 
is 0.0232, with a standard deviation of 0.0094. 
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Figure 4: Calibrated Ho-Lee o parameter over time. 

For the Hull-White model the same procedure is carried out. In the 
Figure [5] (left) we observe the daily values of a and in Figure [5] (right) the 
values of a. The mean value for a is 0.0693 with a standard deviation of 
0.0257. The corresponding mean and standard deviation for a are 0.0177 
and 0.0079 respectively. 

The calibration of the Ho-Lee model corresponding to Wednesday 25th 
of June of 2014 gives a = 0.3071, and is shown in Figure [6] (left). In Figure 
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Figure 5: Calibrated Hull-White parameters over time. On the left a, on 
the right a. 

[6] (right) we see the calibration for the same day corresponding to the Hull- 
White model. The parameters in this case are a = 0.0813 and a = 0.0215. 
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Figure 6: Bond price curve adjusted to market data. On the left Ho-Lee, on 
the right Hull-White. In both cases the dots represent the market prices. 

Both in the Ho-Lee and Hull-White models, the daily calibration allow 
us to construct term structure of interest rates for the corresponding day. 
The result is shown in Figure [7} 

7 Conclusions 

The objective of the present work is to present an arbitrage-free consistent 
model to price the Uruguayan debt nominated in USD. A second purpose is 
to provide an instrument capable of pricing derivatives on bonds, as interest 
rate swaps that are beginning to be used in the Uruguayan market. Based on 
data from the Uruguayan market for coupon bonds we adjust four different 
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Figure 7: Term structure in Uruguayan Market via Hull-White Model. 


gaussian models. 

This information is of valuable interest to financial practitioners and 
policy makers alike. Policy makers monitor expectations of future monetary 
policy to gauge the effectiveness of their strategy. For practitioners, the 
availability of accurate interest rate forecasts can be the key to a success¬ 
ful trading strategy. We hope that this modelling exercises would enrich 
our understanding of market expectations and improve the characteristic 
behaviour of the term structure of interest rate. 

First we adjust the Vasicek and the G2++ models for interest rates, de¬ 
parting from data from one bond and two bonds respectively. We follow the 
methodology proposed by Chen and Scott [3] that departs from historical 
data to estimate the parameters of the model through maximum likelihood. 
The second model is more flexible than the first one, as it involves more 
parameters, and also adjust the initial interest rate curve avoiding the possi¬ 
bility of arbitrage. Nevertheless, as we show, when high values of correlation 
are used (and this is our case), the model can give arbitrage situations. 

Second, we model the debt through the HJM models. We choose two 
different situations, assuming first that the volatility is constant, giving raise 
to the Ho-Lee model, the simplest HJM model, and second we assume a 
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exponentially decreasing volatility, giving rise to the two parameter Hull- 
White model. In this case the adjustment seems to be better, as no arbitrage 
possibilities appear (the HJM model uses the initial curve and also works 
under a risk-neutral measure), and the different values for different days in 
both cases show certain stability. We conclude that HJM models fit better 
to the data, and are capable of describing the whole market structure with 
information of a series of then more traded bonds. 

Further work includes the consideration of more sophisticated HJM mod¬ 
els, some test to verify the adjustment of issued bond prices that are not 
used in the calibration procedures, and the computation of derivatives prices 
in bonds. 
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